Conductance of ID quantum wires with anomalous electron-wavefunction localization 
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We study the statistics of the conductance g through one-dimensional disordered systems where 
electron wavefunctions decay spatially as ~ exp(— \r a ) for < a < 1 , being a constant. 
In contrast to the conventional Anderson localization where ~ exp(— Ar) and the conductance 
statistics is determined by a single parameter: the mean free path, here we show that when the 
wave function is anomalously localized (a < 1) the full statistics of the conductance is determined 
by the average (In g) and the power a. Our theoretical predictions are verified numerically by using a 
random hopping tight-binding model at zero energy, where due to the presence of chiral symmetry in 
the lattice there exists anomalous localization; this case corresponds to the particular value a = 1/2. 
To test our theory for other values of a, we introduce a statistical model for the random hopping 
in the tight binding Hamiltonian. 

PACS numbers: 72.10.-d, 72.15.Rn, 73.21.Hb 



I. INTRODUCTION 

The phenomena of electron wavefunction localization- 
Anderson localization-in a disordered media has brought 
the attention of physicists for several decades*^ Nowa- 
days signatures of localization have been found in dif- 
ferent physical systems. For instance, experiments with 
light, acoustic waves, microwaves, and cold atoms have 
reported evidence of localization^^ 

In the standard Anderson localization problem, elec- 
tron wave functions are localized exponentially in space: 



\ip\ ~ exp(-Ar), 



(1) 



where A can be identified as the inverse of the localiza- 
tion length. For practical purposes, it is more conve- 
nient to define the localization length through measur- 
able transport quantities; for a system of length L, the 
localization length is defined by the exponential decay of 
the dimcnsionlcss conductance g, or transmission. Since 
g oc \ip(L)/ip(0)\ 2 we have that g oc exp (— 2XL). Thus, 
the inverse localization length A is usually estimated by 
the relation 



(-In g) = 2AL, 



(2) 



i.e., the average (lag) is a linear function of L in the 
standard electron localization problem. Within a non- 
interacting electron model, a scaling approach of local- 
ization has successfully described the statistical proper- 
ties of electronic transport Within this approach, it 
has been found that the complete distribution of the di- 
mensionless conductance is determined by a single pa- 
rameter: the inverse localization length 1 ^, given by Eq. 
([2]). In general, one might say that there is a good un- 
derstanding of the statistical properties of the transport 
in the Anderson localization problem in one dimensional 
(ID) and quasi-one dimensional disordered systems. 



On the other hand, anomalous localization of elec- 
tron wave functions has been found in ID disordered 



systems j 



against the general idea that in ID sys- 



tems all the electronic eigenstates are always exponen- 
tially localized. This problem has been much less stud- 
ied than the above standard localization phenomena. For 
instance, a disordered system described by a random hop- 
ping tight binding model was studied in Ref. ITU , where 
it was found that the typical conductance (exp (In g) ) be- 
haves as 



Styp 



exp(-AL 1/2 ). 



(3) 



This unconventional localization of electrons (also named 
delocalization 1 ^) can be explained by the presence 
of a symmetry in the lattice, the so-called chiral 
symmetry , 13 ' 14 which makes the energy spectrum sym- 
metric around zero energy.— The effects of the chiral 
symmetry in a disordered system was studied also within 
a scaling approach to localization! 15 i 16 It was found that 
there is no exponential localization of the conductance 
and the logarithm of g is not self-averaging, while the en- 
semble average (In g) is not proportional to L, as in the 
standard Anderson localization, but to L 1 / 2 , i.e., (In g) oc 
L 1 / 2 . A similar delocalization has been found in disor- 
dered superconducting wires,— ~— where the Bogoliubov- 
de Gennes Hamiltonian has additional symmetries. 21 De- 
localization at zero energy has been also studied using 
tight-binding models of spinless fermions with particle- 
hole symmetric disorder— and in ID systems in the con- 
text of phase transitions in random XY spin chains^ 3 - 
which is mapped onto the so-called random mass Dirac 
model; within this model, it was also foun d 24 ' 25 that 
(—lag) oc L 1 / 2 . In addition, statistical properties of the 
conductance in 2D systems under the presence of chiral 
symmetry has been studied in Ref. H^. 

In the present paper we show that the complete distri- 
bution of the conductance for anomalous transport (non- 
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standard exponential localization) can be determined by 
the value of the average (In g) and the power a of its 
dependence on length L, i.e., (In g) oc L a . Thus, within 
a model of noninteracting electrons, the microscopic de- 
tails of the systems (Hamiltonian) do not enter into the 
description of the statistical properties of the transport, 
in this sense, the description is universal. Our theoret- 
ical model is based on a previous study of the conduc- 
tance statistics of ID disordered quantum wires where 
the random configuration of potential scatterers along 
the wire follows a distribution with a long tail (Levy- 
type distribution) j2! However, in that paper, the anal- 
ysis of the transport was restricted to disordered wires 
where information on the Levy-type distribution was ex- 
plicitly introduced into the disorder configuration of the 
scatterers. Here, we do not need a Levy- type disorder 
configuration but a mechanism to produce anomalous lo- 
calization of the electron wave function, within a single- 
electron model, e.g. the chiral symm etry . Thus, as we 
show in this work, the results in Ref. [27| can be applied 
in general to disordered systems where electron anoma- 
lous localization is present. This larger scope of such 
statistical analysis was overlooked in Ref. 

The remainder of this paper is as follows, after present- 
ing a brief review of the results for wires with Levy-type 
disorder, we introduce the random hopping tight bind- 
ing model where at zero energy anomalous localization 
is present. The numerical results of this model will be 
compared with our theoretical predictions; in particular, 
we arc interested in the conductance distribution. The 
numerical results from the random hopping tight binding 
model at zero energy corresponds to a special case of our 
theory (a = 1/2). To go further and verify our results in 
a more general way, we introduce a statistical model for 
the random hopping which allows to study different de- 
grees of localization characterized by the value of a. We 
finally summarize our results and give some conclusions 
in the last part of the paper. 



processes have been studied in connection to transport 
problems). It was found in [27| that the full statistics of 
the conductance is determined by the average (In g) and 
the exponent a of power-law tail in the macroscopic limit 
(L c x / a ). In particular, it was shown that the com- 
plete distribution of conductances P^(g), with £ = (ln g), 
is given by 

/•OO 

P${9) = / Ps( a .t. z )(g)q a ,i(z)dz, (4) 

J 

for a < 1, where q a ^ c is the probability density function 
of the Levy-type distribution supported in the positive 
semiaxis, s(a, £, z) — £/ (2z a I a ), I a = 1/2 j Q z~ a q a ,idz, 
and 

s 2 e 4 / we 4s 

where yo = arcosh(2/g — 1). Also, it was shown that the 
average of the logarithm of the conductance depends on 
L as 

(-lng)cxL Q , (6) 

for < a < 1, while for values 1 < a < 2 the linear be- 
havior ((— In g) oc L) is recovered. From the same model 
one can also find that conductance average behaves as 

(g) « L~ a , (7) 

for < a < 1, in contrast to the exponentially depen- 
dence with L in the standard localized regime. The most 
interesting effects of anomalous localization are seen for 
values < a < 1, so we concentrate in this region, al- 
though the case 1 < a < 2 can be analyzed within the 
same theoretical framework. 

III. ANOMALOUS LOCALIZATION: a = 1/2 



II. THEORETICAL MODEL 

As we have mentioned, our theoretical model of this 
work is based on a study of coherent transport in the 
presence of Levy-type disorder^ We briefly mention 
that Levy-type random processes are described by a den- 
sity probability q UiC (x) with a long tail: for large x, 
Qa,c(x) ~ c/x 1+a with < a < 2 and c being a constant. 
These kind of distributions are also known by mathe- 
maticians as a-stable distributions i22r— Notice that first 
and second moments diverge for < a < 1. Motivated 
by the realization of experimentally controlled Levy pro- 
cesses in the so-called Levy glasses^ in Ref. [27[ a model 
was developed to describe the statistical properties of the 
conductance through a ID quantum wire where electrons 
suffer multiple scattering due to scatterers placed along 
the wire in a random way accordingly to a Levy-type 
distribution (see f33U37j for other examples where Levy 



Next we consider the tight binding model with nearest 
neighbor random hopping, at zero energy, described by 
the Hamiltonian 

H = ^ tn(clc n+1 + c,\ +1 c„), (8) 

n 

where and c n are creation and annihilation operators 
for spinless fermions, and t n (> 0) are the random hop- 
ping elements sampled from a distribution of the form 
Pit) — l/wt,exp(—w/2) < t < exp(w/2), where w 
denotes the strength of the disorder. This is the so- 
called logarithmic off-diagonal disorder^ As we have 
mentioned, the model described by Eq. (j8|) has been 
found to present unconventional localized states at zero 
energy ^ ~ 14 i 22 ' 24 i 25 whereas for nonzero energy standard 
localized states are present. To illustrate this fact, we 
have calculated the conductance within the Landaucr- 
Buttiker approach. In Fig. Q] we show the ensemble av- 
erage (In g) as a function of the length of the system (in 
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FIG. 1: ( — In g) as a function of the length L at energy E = 
for strength of disorder w = 2.5. The solid line is obtained 
by fitting the data (dots) according to Eq. (|6]) with a = 1/2. 
Upper inset :{g) as a function of the length L (for the same 
parameters in the main frame). The solid line is fitted to the 
numerical data assuming that (g) oc L _1//2 . A good agree- 
ment is seen. Lower inset: (~\ng) for a linear chain with 
off-diagonal logarithmic disorder as a function L at energy 
E — 0.1 and strength of disorder w — 2.5 (50000 realiza- 
tions). As expected, a linear behavior is observed indicating 
Anderson localization. 




-10 -8 -6 -4 -2 o 



FIG. 2: The distribution of In g for a system of length L = 400 
with offdiagonal logarithmic disorder, at energy E = and 
strength of disorder w — 0.35 (50000 realizations). From the 
numerical data {—\ng) = 2.1. Using this information the 
theoretical distribution (solid line) is calculated with a = 1/2, 
Eqs. ([4]) and (J5]). Inset: P(g) for the same case as in the main 
frame. The coexistence of insulating and ballistic regimes are 
manifested by the presence of two peaks at g — and 1. As 
we can see, the theory (solid line) gives correctly the trend of 
the numerical results (histograms). 



units of the lattice constant) at zero and nonzero ener- 
gies. As we can observe (— In g) oc L 1 / 2 at zero energy 
(main frame) , while a linear dependence on L is obtained 
at finite energy (lower inset), restoring the standard An- 
derson localization. Additionally, in the upper inset of 
Fig. Q] we show the average of the conductance (g) at 
zero energy, which depends on the length as L^ 1 / 2 , as 
given by Eq. J7]). 

We now show that the complete distribution of con- 
ductance is described by Eq. As we have claimed, in 
order to compare the theoretical and numerical results, 
we only need the information of the value (In g) and its 
power dependence on L, which are taken from the nu- 
merical simulation; thus, there is no free parameters in 
our theory. In Figs. [2] and |3] we show the distribution of 
the conductance obtained from the numerical simulations 
(histograms) for two different strengths of disorder and 
the corresponding theoretical distributions (solid lines) 
accordingly to Eq. dU). Note that we plot P(lng) in 
the main frames, instead of P{g), since for very insulat- 
ing cases the details of the distributions are better seen 
in this way. For the smaller case of strength disorder 
(w = 0.35) in Fig. [2] we have included P(g) in a inset. 
Here we can observe two peaks at g = and g = 1, 
which is due to the existence of strong sample-to-sample 
conductance fluctuations, i.e., in our ensemble a consid- 
erable amount of samples behaves like insultors {g << 1), 
whereas another important amount of them behaves as 
ballistic samples (g ~ 1). This behavior is very robust 
in the sense that if we increase the length of the system 
or the disorder degree the peak at g = 1 survives. This 
is not seen in the conventional ID electron localization 
problem. In Fig. [3] we increase the strength of disorder 
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FIG. 3: The distribution of lng for strength of disorder 
w = 1.2. (—big) — 9.7 corresponding to a more insulating 
case than the previous one (Fig. [2]), while the power-law 
dependence on L remains 1/2 (Fig. [1]). A good agreement is 
seen between the numerical histogram and the corresponding 
theoretical distribution (solid line). 



to w = 1.2. Thus, for both strengths of disorder, Figs. [5] 
and [3] show that our theory gives correctly the trend of 
the numerical distribution. We might see a small differ- 
ence between numerics and theory in the inset of Fig. [2] 
at g « 1, but we would like to remark that there is not 
free parameter in our theory. Therefore our model with 
a = 1/2 describes correctly the statistics of the conduc- 
tance when anomalous localization of the wave function 
is of the form \ip\ ~ exp (— XL 1 / 2 ). However, this is a spe- 
cial case for our model. We would like to explore different 
exponential power decays a of the wave function. 
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FIG. 4: (— In g) and (g) (inset) as a function of the variable 
T for a = 1/3 in the statistical model for the tight bind- 
ing Hamiltonian (see text). 50000 realizations are considered 
and energy E — 0.1. The solid lines are obtained by fitting 
the power dependence: T 1 ' 3 and T _1//3 for (— Ing) and (g), 
respectively, as predicted by the theoretical model. 



IV. ANOMALOUS LOCALIZATION: 
ARBITRARY a 

In order to investigate different anomalous-localization 
degrees of the wave function, we introduce a statistical 
model for the nearest-neighbor random hopping model, 
Eq. ([5]). In fact, what we need is a model that induce 
large fluctuations of the conductance. A way to introduce 
such a large fluctuations is to consider the hopping t n 
as a random variable that follows a distribution with a 
long tail, i.e., a Levy-type distribution, and keeping fixed 
the total sum of the hopping elements: T — t n . By 
varying the value of T we can change the degree of the 
localization of the disordered samples. We have verified 
numerically that T acts similarly to the length L in the 
Levy-type configurational disorder used in [23. However, 
the tight binding model is more appropriate for numerical 
simulations. The study is carried out at non-zero energies 
in order to get rid of the effects of chiral symmetry. 

With the above statistical model for the random hop- 
ping tight binding Hamiltonian, we calculate the statis- 
tics of the conductance. The data are collected over an 
ensemble of 50000 realizations of disorder. In Figs. [Hand 
[5] we show first the results for the average (— In g) and 
(g) (insets) as a function of T where the random hopping 
elements are generated from two different Levy-type dis- 
tributions with tail decay exponents a =1/3 and 3/4. 
We can see that indeed (— In g) cx T a and (g) cx T~ a , for 
both values of a. Having in mind that T plays a similar 
role as L in our configurational disorder model^l we ex- 
pect that the wave function is anomalously localized as 
|V>| ~ exp {-XL 1 / 3 ) and ~ exp (-XL 3 / 4 ), for a = 1/3 
and 3/4, respectively. 

We now show that the distribution of the conductance 
is described by Eq. (HJ). For a — 3/4 and two different 
values of T, in Figs. [H] and [7] we compare the numerical 
simulations (histograms) and the corresponding theoret- 
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FIG. 5: Numerical data (dots) of (— lng) and (g) (inset) from 
the tight binding model at energy E = 0.1 and a — 3/4. 
The solid lines are fitted assuming that (—lng) oc T 3 / 4 and 
(g) oc T~ 3//4 , in agreement with the model, Eqs. ([6| and lf?|l. 




FIG. 6: The numerical distribution P(\ng) (histograms) for 
a = 3/4 with E = 0.1, T = 35, and {-lng) = 2.0 . Inset: 
P(g) for the same parameters of the main frame. The solid 
line is obtained accordingly to Eq. . A good agreement 
between numerical and theoretical (solid line) results is seen. 

ical results (solid line). The case in Fig. [6] is less insu- 
lating than that one in Fig. [7J so we plot in an inset 
the distribution P(g). For the more insulating case (Fig. 
[Jj , we can observe a nonconventional shape of the distri- 
bution P(hig). We mean by nonconventional shape the 
non Gaussian shape of the distribution; we recall that 
for the standard Anderson localization it is expected a 
log-normal distribution in the insulating regime. Thus, 
from both Figs. [S] and [7] we can see that the trend of the 
numerical distributions are well described by our theory. 
Finally in Fig. [8] we show the distribution P(\ng) for 
a = 1/3. Here we also note the nonconventional shape 
of the distribution which is a consequence of the anoma- 
lously large conductance fluctuations. 

V. CONCLUSIONS 

To conclude, in this work we have shown that the com- 
plete statistics of the conductance of an ID disordered 
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FIG. 7: The distribution of hip for a = 3/4, T = 250 and 
E = 0.1. For this case {— \ng) = 9.3. We can see that the 
theoretical result (solid line) describes correctly the numerical 
distribution. 
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FIG. 8: The distribution of In g for a = 1/3, T = 100 and 
energy E — 0.1. {— In g) = 9.4 for this case. Comparison with 
the corresponding theoretical distributions is shown. A good 
agreement between theory and numerics is seen. 



system, when electron wave functions are anomalously 
localized ( ip ~ cxp (— Xr a ), < a < 1), is determined 
by the exponent a and the average (In g). In contrast, 
in the standard Anderson localization, the knowledge of 
(Ing) is enough to describe the statistical properties of 
the conductance. We have verified our results for differ- 
ent values of a. For the particular case of a = 1/2, we 
have used a random hopping tight binding Hamiltonian 
at zero energy to verify our predictions since it is well 
known that nonexponential localization in this model is 
present due to the existence of chiral symmetry on the 
lattice. In order to study other degrees of anomalous 
localization (different values of a) we have introduced 
a statistical model for the hopping in a tight binding 
Hamiltonian that promote the presence of large fluctua- 
tions of the conductance. We remark that our theoretical 
model do not make any reference to a specific Hamilto- 
nian system and there is no free parameter; the infor- 
mation needed in our theoretical model (a and (lng)) is 
extracted from the numerical simulation. On the other 
hand, we have restricted our study to ID systems (one 
channel) , we think an extension to multichannel systems 
is of interest since other regimes of transport, e.g. the 
diffusive regime, can be analyzed. Finally, the conduc- 
tance statistics in the conventional Anderson localization 
problem has been extensively studied, we hope this work 
helps to the understanding of a much less studied topic in 
quantum transport: the statistical properties of the con- 
ductance when electron wave functions are anomalously 
localized. 
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